ENSLAPP-A-643/97 



The Production of Anti-Matter in our Galaxy 

Pascal Chardonnet'^''', Jean OrlofP^ and Pierre Salati'^''^ 



a) Laboratoire de Physique Theorique ENSLAPP, BPllO, F-l^O^l 
Annecy-le-Vieux Cedex, France. 

b) Universite de Savoie, BPIIO4 73011 Chambery Cedex, France. 



Abstract 

The discovery of a single anti-helium nucleus in the cosmic ray 
flux would definitely point toward the existence of stars and even of 
entire galaxies made of anti- matter. The presence of anti- nuclei in 
cosmic rays has actually profound implications on the fundamental 
question of the baryon asymmetry of the universe. It is therefore 
crucial to determine the amount of anti-matter which our own galaxy 
already produces through the spallation of high-energy protons on the 
interstellar gas of the galactic disk. We have used here a coalescence 
model to assess the amount of anti-deuterium and anti-helium ^He 
present in cosmic rays together with anti-protons. The propagation 
of cosmic rays in the galaxy is described through a two-zone diffusion 
model which correctly describes the observed abundances. We find 
that the D/p ratio exceeds 10~^ above a momentum per anti-nucleon 
of ~ 4 GeV/c. Would the universe be purely made of matter, the AMS 
collaboration should be able to detect a few anti-deuterons during the 
space station stage of the experiment. However, the ^He/p abundance 
does not exceed ~ 4 x 10~^^. Heavier anti- nuclei are even further 
suppressed. 



1 Introduction. 



The amount of anti-matter in cosmic rays is about to be measured with un- 
equalled accuracy by the space shuttle borne spectrometer of the AMS collab- 
oration One of the most exciting goals of the experiment is the possible 
detection of anti-nuclei in the cosmic radiation. It is generally believed that 
the observation of a single anti-helium or anti-carbon would undoubtedly 
signal the presence of stars made of anti-matter. Such a discovery would be 
of paramount importance as regards the existence of a baryon symmetry in 
the universe and has therefore strong cosmological implications. 

That is why it is crucial to ascertain that cosmic rays do not already 
contain detectable traces of anti-nuclei which could have been directly man- 
ufactured in our galaxy. We know for instance that a fraction p/p ~ 2 x 10~^ 
of anti-protons is produced by the spallation of cosmic ray protons on the 
interstellar gas of the galactic disk. In this letter, we compute the abundance 
of anti-deuterium D and anti-helium ^He produced through the interaction 
of high-energy protons with the interstellar hydrogen. 

That calculation requires two ingredients. First, we need to evaluate 
the production cross section of anti-nuclei during the interaction of a high- 
energy proton with a proton at rest. Anti-deuterons are formed by the fusion 
of an anti-proton and anti-neutron pair. We have used here a factorization 
scheme together with a coalescence model [0 which we discuss in section ^. 
Our estimates of the probability for anti-deuterium production are in fairly 
good agreement with the accelerator data collected at Serpukhov and at the 
ISR 0, ^. The anti- helium ^He is predominantly formed through the 
production of anti-tritium T which subsequently decays into ^He. Then, in 
section ^, we recall the salient features of the two- zone diffusion model that 
takes care of the propagation of anti-nuclei in the galaxy from the production 
regions to the earth. Finally, our estimates of the cosmic ray abundances D/p 
and ^He/p are discussed in section ^. We find that the D/p ratio exceeds 
10~^ above a momentum per anti-nucleon of ~ 4 GeV/c. Would the universe 
be purely made of matter, the AMS collaboration should still be able to 
detect a few anti-deuterons during the space station stage of the experiment. 
However, the ^He/p abundance does not exceed ~ 4 x 10~^^. Heavier anti- 
nuclei are even further suppressed. 
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2 Factorization and coalescence. 



Our goal at this point is to establish the Lorentz- invariant cross section for the 
production of the anti-nucleus x during the interaction between two protons. 
Unless stated otherwise, the subsequent reasoning will be performed in the 
center of mass frame of the proton-proton collision with total available energy 
a/s and total cross section cf^^I^- The number dM^ of species x created during 

— * 

a single interaction and whose momenta are is related to the corresponding 
differential probability of production by 

dJ\f^ = J^^{^s,k^)d% . (1) 

This probability JF^ may be expressed in turn as a function of the Lorentz- 
invariant differential cross section 

d^a 

^(v^, K) = ^x(v^> K) <-p , (2) 

which we want to evaluate. The anti-nucleus energy cannot exceed the 
upper bound 

3-^1 + ml 



where Mx denotes the minimum mass of the nucleon system (X) that is 
created with the anti-nucleus. As a result of the conservation of the baryon 
number, each anti-nucleon is actually forced to be produced together with 
an additional nucleon. When an anti-proton is created, for instance, three 
protons remain in the final state and Mx = 3mp. As regards the anti- 
deuterium production, Mx — 4mp while it is 5mp for anti-tritium. 

The invariant cross section for the production of anti-protons is experi- 
mentally well known. It is fairly well fitted by the Tan and Ng's parame- 
terization 1^ which we have used throughout this analysis. Assuming that 
the invariance of isospin holds, the anti-neutron production cross section is 
equal to its anti-proton counterpart. The calculation of the probability for 
the formation of an anti-nucleus can now be performed in two steps. We first 
need to estimate the probability for the creation of a group of anti-nucleons. 
Then, those anti-nucleons fuse together to yield an anti-nucleus. Let us con- 
centrate on the case of anti-deuterons which requires the formation of both an 
anti-proton and an anti-neutron. We have assumed that factorization holds 
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at this stage. This means that the production of two anti-nucleons is pro- 
portional to the square of the production of one of them, a hypothesis which 
is reasonably well established at high energies. However, at lower energies, 
factorization has to break down, if only to respect the kinematic constraints 
that thresholds are different. We propose to take this into account by assum- 
ing in addition that the center of mass energy available for the production of 
the second anti-nucleon is reduced by twice the energy carried away by the 
first anti-nucleon 

^p,n{Vs,kp,kri) = ^J^p{\/s,kp)J='fi{y/s-2Ep,kfi) + (kp^k^ . (4) 

This probability has the merit of vanishing at the physical threshold. How- 
ever, it introduces the arbitrary Ansatz that the first pair is produced back 
to back in the center of mass frame of the proton-proton collision. To test 
how crucial this assumption really is in the final results, we also give plots 
for a second possible Ansatz, where the nucleon associated to the first anti- 
nucleon is produced at rest. This simply amounts to replace ^/s — 2Ep by 
^/s — rrip — Ep in the previous formula. As may be already guessed from 
figure both Ansaetze give similar results. In the plateau regime where 
the anti-nuclei abundances are fairly independent of the energy, the relative 
difference between the two assumptions does not exceed 4%. The probabil- 
ity for the production of the p — li pair is related to the corresponding cross 
section by 

^p; ^n) — p,ii("\/^5 ^n) • (5) 



V'^p-p/ d^kpd^kn 

Although ^'s are not Lorentz- invariant, this relation allows to check the 
Lorentz invariance of the factorization Ansatz (^). 

Once the anti-proton and the anti-neutron are formed, they combine to- 
gether to give an anti-deuteron with probability 

•^d(v^) ^d) '^'^^D ~ J d^kp d^ki^ C(kp, ki^) J-p^fi (^y/s, kp, kn^ . (6) 
The summation is performed on those anti-nucleon configurations for which 

kp -\- kf^ = kf) . (7) 

The coalescence function C{kp, fcn) describes the probability for a p — ii pair 
to yield by fusion an anti-deuteron. That function depends actually on the 
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difference kp — = 2A between the anti-nucleon momenta so that relation 
(H) may be expressed as 

^d(v^, kl) = I d^A C(A) ^p,5 ["v^, kp = ^ + A,h = Y-^ 

(8) 

Notice that the formation of the anti-nucleon pair is implicitly assumed to 
be independent of the later coalescence process where the anti-nucleons melt 
together to form the anti-deuteron 0. This is well justified by the large 
difference between the binding and pair creation energies. As a matter of 
fact, an energy of ~ 3.7 GeV is required to form an anti-deuteron whereas 
the binding energy of the latter is i? ~ 2.2 MeV. The coalescence function is 
therefore strongly peaked around A = and expression (|^) simplifies into 

J't,iV~s,kl) ^ {/ rf3AC(A)} ^p,5 (^V~s,h = = , (9) 

where the probability for the formation of the p — li pair has been factored 
out. From the term in brackets, we build a Lorentz-invariant quantity which 
we evaluate this time in the rest frame of the anti-deuteron 



d^ACiA) ^ l-^ i-nPLA . (10) 
EpEf, \mpmnj V3 / 

In that frame, the anti-nucleons merge together if the momentum of the cor- 
responding two-body reduced system is less than some critical value Pcoai- 
That coalescence momentum is the only free parameter of our factorization 
and coalescence model. The invariant cross section for anti-deuterium pro- 
duction in proton-proton collisions may finally be expressed as 

(11) 

2<ip d% -^'^ - d^k, ^ ' " + ^ ^ ■ 

with kp = fcfi = kf)/2. The dependence on Pcoai comes only from the multi- 
plicative Pcoai prefactor. 




Theoretical values for Pcoai range from yrrip i? ~ 46 MeV, naively derived 
from the anti-deuteron binding energy, up to 180 Mev as would follow from 
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a Hulthen parametrization of the deuterium wave function [0 . We therefore 
expect Pcoai to he somewhere in the range between 50 and 200 MeV. Inside 
this range, since factorization might also involve an unknown coefficient that 
could be reabsorbed into Pcoai, we will rather follow a phenomenological 
approach. We will determine Pcoai directly from the known experimental 
results. These have been summarized in terms of Pcoai on figure ^ Points 1 
and 2 are from an experiment carried out at Serpukhov [^], for a center of 
mass (CM) energy ^/s = 11.5 GeV, at practically vanishing longitudinal CM 
momentum and at large transverse momentum with p± = 1.15 and 1.5 
GeV respectively. Points 3, 4 and 5 have been obtained at the ISR 0, with 
^ = 53 GeV, smaU p± (0.16, 0.21 and 0.30 GeV) and larger pl (4.8, 5.7 and 
8 GeV). Finally point 6 summarizes about 8 ISR data points at vanishing 
Pl and with p±_ ranging from 0.2 to 1 GeV. Given the crudeness of our model 
and the wide range of kinematic regimes explored experimentally, it is quite 
comforting that all these data are compatible within 2 standard deviations 
with our predictions from our single parameter model. Notice in that respect 
that most of the error bars simply ignore systematic errors. The first five data 
points are all fitted with a coalescence momentum of order 60 MeV. Because 
the Serpukhov data correspond to a low center of mass energy, a kinematic 
regime where anti-deuteron production is astrophysically predominant, we 
have decided to fix the value of Pcoai at 58 MeV. Increasing this value to 75 
MeV would simply double the production of anti-deuterium. 

The anti-helium ^He nucleus may be formed directly by the fusion of 
two anti-protons and an anti-neutron. A second possibility is the synthesis 
of anti-tritium T which subsequently decays into ^He with a half-lifetime 
of ~ 12.3 years. This second mechanism is dominant. It does not suffer 
from the electromagnetic suppression due to the Coulomb repulsion between 
the two anti-protons. We have therefore ignored here the direct ^He fusion 
and have concentrated instead on the production of T. The corresponding 
production cross section may be obtained from the direct generalization of 
the anti-deuteron calculations. This yields 

^ iUpiy (^] X (12) 





' rrif ^ 


d^kf ~ ' 









3--coaiy ^^tot^^ 



(v^ - 2Ep, Q Ep^ (yi - 4Ep, h] 



d^ kp kp 



The momenta /ci, k2 and /cs of the three anti-nucleons are all equal to /cf/S 
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so that their energy is also the same. Unfortunately, there are no data 
available to calibrate the coalescence factor in this case. The best we can 
do is to use the value Pcoai = 58 MeV extrapolated from the anti-deuterium 
data. We have checked that relation (|1^) does not violate the experimental 
upper bounds f^. 

3 The diffusion of cosmic rays in the galaxy. 

Cosmic rays are produced when supernovae shocks accelerate the interstel- 
lar material of the galactic plane. These high-energy nuclei propagate in the 
erratic magnetic fields of the disk where they interact on the gas to create sec- 
ondary cosmic rays. The amount of secondary light elements such as lithium, 
beryllium and boron (Li-Be-B) is well explained by the spallation of primary 
carbon, oxygen and nitrogen nuclei. The latter spend ~ 10 million years 
(My) in the galactic plane where they cross a column density of ~ 15 g cm~^. 
Quite exciting is the measurement of the isotopic ratio between the unstable 
^°Be and its stable partner ^Be. The former nucleus has a half-lifetime of 
1.6 My and plays the role of a chronometer. A low value is observed and 
indicates that cosmic rays spend actually 100 My in the galaxy before escap- 
ing in outer space. Particles are confined 90% of the time in extended layers 
above and beneath the matter ridge where they just diffuse without inter- 
acting much with the scarce interstellar medium. Therefore, as regards the 
propagation of cosmic rays, our galaxy may be reasonably well modelled with 
two main regions. First, cosmic rays are accelerated within a thin disk with 
radius < r < R = 20 kpc and thickness \z\ < h = 100 pc, where they 
diffuse and interact on atomic and molecular hydrogen. This gaseous plane is 
sandwiched by extended domains containing irregular magnetic fields, with 
same radial extension and \z\ < L = 3 kpc. These thick layers play the 
role of confinement reservoirs. As a matter of fact, the presence of magnetic 
fields far above the galactic plane is now firmly established by synchrotron 
radiation as mentioned by Badwar [Q. On average, the intensity of the mag- 
netic field decreases from ~ 5 fiG in the plane down to 1 /iG at a height of 5 
kpc. Note that magnetic fields have also been detected in other galaxies 0. 
Therefore, the transport of high-energy particles crucially depends on their 
diffusion across the chaotic magnetic fields of the galaxy. We will assume 
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here an isotropic diffusion with coefficient given by the empirical value 



K = 6x10^3^2^-1 + , (13) 

where TZ = p/Z is the particle rigidity. The diffusion coefficient K is constant 
at low energies, but increases with the rigidity beyond a critical value of 3 
GV. 



We follow here the analysis of Webber, Lee and Gupta |]T0[- The au- 
thors showed that the two-zone diffusion model presented above is in good 
agreement with the observed primary and secondary nuclei abundances. As- 
suming that steady state holds, the distribution of cosmic ray protons may 
be derived from the diffusion equation 

- = (14) 



dt 



where np{Ep, r, z) is the density of protons with energy Ep at location (r, z). 
In the right-hand side of relation ([T^, the first term describes the diffusion 
of particles. The second term accounts for the destruction of protons by 
their interactions with the interstellar medium of the galactic plane. That 
destruction rate is shown to be very small but has not been neglected in 
what follows. The disk is only 200 pc thick, hence our approximation of an 
infinite thinness and the term 2h6{z). The collision rate of protons with the 
interstellar hydrogen is given by 

Tp = al% Vp riH , (15) 

where (J^lp ~ 44 mbarns is the total proton-proton interaction cross section, 
Vp denotes the velocity and nn = 1 cm"^ is the average hydrogen density in 
the thin matter disk. The last term in relation (|T^) deals with the sources 
of high-energy protons. It matches the distribution of supernovae remnants 



and pulsars as measured by Lyne, Manchester and Taylor [|TT 

Qp{Ep,r) = QliEp) p'-' e-'P , (16) 

with p = r/R. At the edge of the domain where cosmic rays are confined, 
the particles escape freely and diffusion becomes inefficient. Thus the density 
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vanishes at the boundaries of the axi-symmetric domain where cosmic rays 
propagate. This provides the boundary conditions for solving the diffusion 
equation (|I^). Following |T0[, the density Up is expanded as a series of the 



Bessel functions of zeroth order Jq {(ip) where Q is the ith zero of Jo and 
where p = r/R. At fixed proton energy Ep and corresponding rigidity Tip, 
the density may be expressed as 

oo 

np{r,z) = Y.Np^{z)Jo{C^p) ■ (17) 

i=l 

The diffusion equation ([1^) may be Bessel transformed, yielding 

K - ^NpA - 2h5{z)TpNp, + 5{z)q, = . (18) 

The coefficients qi = 2h Qi are the Bessel transforms of the source distribution 
Qp{r) in the galactic plane 



[Jl^^C^)}[[Qp{p)Jo{C^p)pdp] [j\Mpdp 



(19) 



where Qp°* denotes the total production rate of cosmic ray protons with 
energy Ep in the entire galactic ridge. We have assumed here that the energy 
spectrum of the sources has the same shape all over the disk. Therefore, the 
energy dependence of the coefficients may be factorized out in the term 
Q^°^{Ep) and does not depend on i. The functions Npi are derived from the 
requirement that they vanish at the boundary z = L and from the relation 

2K 



dz 



= 2hrpNp,{0) - q, . (20) 

z = 

Integrating equation (|THp across the disk and remembering that the A^pj's 

are odd functions of the height z leads to the last expression. The cosmic 
ray proton density may finally be expressed as 

A^P.(^) = -j^H^) ■ (21) 

The vertical distribution J-'i{z) is given by 

sinh I — (i 

Hz) = ^\r.. ^ ' > (22) 



sinh < 


>- 


A)] 


sinh < 




\ 
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where the parameters are equal to K^ijR. The coefficients A^i stand for 

A^i = 2hTp + Kaicothl^^lj , (23) 

and depend on the proton energy through the diffusion coefficient K(Jlp) 
of relation ([I3D and, in a lesser extent, through the collision rate Fp. These 
coefficients are actually dominated by the diffusion contribution. In the case 
of cosmic ray protons, spallation does not play an important role in deter- 
mining their galactic distribution. It affects heavier nuclei for which the 
interaction cross section on hydrogen is significantly larger. 

The density of the anti-nucleus species x may be readily inferred from 
a similar reasoning. It satisfies the diffusion equation 

KAn^ - 2hS{z)T^n^ + 2h5{z)Q^{r) = . (24) 

The source term Qx{f) obtains now from the convolution of the proton energy 
spectrum rip with the cross section for production of an anti-nucleus x ^ 
proton-proton interaction 

Q,(i?„r) = y^^^^jpnp(i?p,r,z = 0) |^(i?p ^ dE, . (25) 

The energy threshold above which the cosmic ray proton may produce an 
anti-nucleus x depends on the atomic number of the latter. For anti-protons, 
E"™™ = 7mp while it is 17 rrip for anti-deuterons and 31 mp for anti-tritons. 
The differential energy distribution dM-^jdE-^ corresponds to the anti-nucleus 
X with energy E^ that is produced when a proton with energy _Ep collides on 
an hydrogen atom at rest. It is related to the Lorentz- invariant production 
cross section estimated in section ^ 

where corresponds to the process p + p ^ x + X. The momentum of 
the outgoing species x makes an angle Q with respect to the direction of the 
incoming cosmic ray proton. The production term may be expanded in 
terms of its Bessel transforms 5*^. This yields the anti- nucleus distribution 

°° ( 2h S'A 

= E H^)Jo{Q.p) , (27) 
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where 



B^i{E^) = 2hT^ + ir(7^JaiCoth|^L| . (28) 

Relation (^) may be simplified because the Bessel transforms N^i of the 
proton distribution rip have all the same energy behavior. The latter is 
contained in the factor Q*^^ (Ep) / K [Ep) as may be inferred from expression 
([2l|) . An effective anti- nucleus multiplicity may therefore be defined as 



No convolution with the energy spectrum Hp of the cosmic ray protons is 
needed any more. The differential flux of protons with energy Ep is denoted 
here by 

%{Ep) = ^UpVp , (30) 

and is expressed in units of cm~^ s~^ sr~^ GeV~^. In order to derive that 
flux, we could have naively computed the proton density Up from the Bessel 
expansion (0), provided Qp{Ep) is well known. However, as explained above, 
cosmic ray protons are produced with the same energy spectrum all over the 
galactic ridge. Their propagation is also dominated by mere diffusion and 
not spallation on the interstellar gas. Therefore, their energy spectrum does 
not depend on the position in the galaxy. The energy dependence of Up 
may actually be factorized. As a consequence, the flux ratio that appears in 
the effective multiplicity (p9D is constant throughout the galaxy and may be 
directly derived from measurements performed at earth. We borrowed the 



proton spectrum from 



<l>p{Ep) (X , . (31) 

+ (1.5GV/7^p)' 

The distribution of the anti-nucleus species x finally be expressed as 



n^{E^,r,z) = E^.(z)Jo(Op) —^2hTpAff . (32) 

i=i ^p«-°x« 



Note that the anti-nucleus and proton energies are both equal to E^. in the 
last formula. The diffusion coefficient K that appears in the quantities Api 
and B^i is therefore different for the proton and for the anti- nucleus because 
at fixed energy, the corresponding rigidities are different. Finally, the energy 
dependence Qp{Ep) in the proton sources qi disappears in the ratio n-^/up. 



10 



4 Discussion and conclusions. 



The flux of cosmic ray anti-protons relative to the flux of protons is plotted 
(solid line) in figure @ as a function of the momentum of the particles. The 
curve drops sharply below a few GeV. At higher energies, it exhibits a plateau 
where the ratio p/p reaches a value of ~ 2 x 10~^, in good agreement with the 
previous calculation by Gaisser and Schaefer [0]. The D/p and ^He/p cases 
are also presented on the same plot. The energy behavior is fairly the same 
as for the anti-protons. However, the magnitude of the effect is significantly 
suppressed. To fit on the same diagram, the curves have been scaled by a 
bench-mark factor of 10'* for anti-deuterium and of 10^ for anti-helium ^He. 
As explained in section 0, the prescription that sets the energy available to 
the production of any additional anti-nucleon in a proton-proton collision is 
not clear. Two options are presented here where the nucleon associated to a 
final anti-nucleon is either produced back to back with its anti-partner (lower 
curve) or at rest (upper curve). The effect is negligible since the magnitude 
of the various plateaux is not affected. 

We find that the D/p ratio exceeds 10~^ above a momentum per anti- 
nucleon of ~ 4 GeV/c. It reaches a maximum of 6 x 10~^ for a momentum 
of ~ 20 GeV/c. Note that we have only computed the anti-deuterium yield 
resulting from the spallations of cosmic ray protons on the interstellar hydro- 
gen. Anti-deuterons may alternatively be produced when anti-protons this 
time interact on the galactic gas. We nevertheless think that this effect is 
negligible. On the one hand side, the pp interaction needs to be quite elastic 
to allow for the survival of the incoming anti-proton. On the other hand side, 
if the collision is not inelastic, the additional anti-neutron is not produced. 
Most of the time, the anti-proton annihilates. The probability that a second 
anti-nucleon is created with same momentum as the incoming anti-proton 
is therefore extremely small. In addition, the anti-proton flux is 4 orders of 
magnitude fainter as for protons. 

During the space station borne stage of the experiment, the AMS collab- 
oration should measure a flux of anti-nucleons with a sensitivity reaching a 
level of ~ 10~^ relative to cosmic ray protons. The corresponding energy per 
nucleon ranges from a GeV up to 20 GeV. We therefore conclude that AMS 
should detect a few cosmic ray anti-deuterons. 

The ^He/p abundance does not exceed ~ 4 x 10~^^. Even allowing for 
a generous error of 2 in our estimates, we conclude that the anti-helium 
nuclei that are manufactured in our galaxy will not be detected by AMS. 
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Note that we have actually computed the interstellar flux. The results of 
figure ^ do not include solar modulation. Since the various curves reach 
their maxima for a momentum per nucleon in excess of 10 GeV, the effect 
of solar modulation will just be a slight shift of a fraction of a GeV towards 
low energies, depending on the specific epoch of the solar cycle. The same 
conclusions still hold. 

As shown in section 0, factorization comes into play for the production 
of several anti-nucleons during the same proton-proton event. The general 
trend is a reduction of the yield by a factor larger than 10^ when the atomic 
number of the anti- nucleus is increased by one unit. As an illustration, it is 
clear from figure ^ that the D/p ratio is ~ 3 x 10^^ while ^He/D ~ 7 x 10^^. 
The anti-nuclei with atomic number A > 3 that are produced through cosmic 
ray spallations in our galaxy are not detectable by AMS. Alternatively, if anti- 
matter is processed in anti-stars, the abundance of its various components 
should be quite similar to the conventional stellar yields. It would follow 
therefore the same relative proportions as in our own interstellar material. 
In any case, the detection of a single anti-helium or anti-carbon by the AMS 
collaboration would be a smoking gun for the presence of large amounts of 
anti-matter in the universe and for the existence of anti-stars and of anti- 
galaxies. 
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Figure 1: This figure displays various experimental constraints on the coa- 
lescence momentum Pcoai, the only free parameter of the model discussed in 
section (^. Points 1 and 2 are from Serpukhov with ^/s= 11.5 GeV while all 
the other data have been collected at the ISR at y/s = 53 GeV. A coalescence 
momentum Pcoai of order 60 MeV provides a reasonable fit of all the points 
but the last one. 



14 



pbar/p 

dbar/p * 1e4 (ansatz 1) 

dbar/p * 1 e4 (ansatz 2) 
heSbar/p * 1 e8 (ansatz 1 ) 
heSbar/p * 1 e8 (ansatz 2) 
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Figure 2: The fluxes of cosmic ray anti-protons and of anti-deuterium and 
anti-tritium nuclei, relative to the proton flux, are presented as a function 
of the momentum per nucleon. To fit on the same diagram, the curves have 
been scaled by a factor of 10'^ for anti-deuterium and of 10^ for anti-helium 
^He. The doubling of curves corresponds to different factorization schemes 
as explained in section |l|. 
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